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Abstract 

The effects of both density stratification and magnetic field expansion on torsional Alfven 
waves in magnetic flux tubes are studied. The frequencies, the period ratio P1/P2 of the 
fundamental and its first-overtone, and eigenfunctions of torsional Alfven modes are obtained. 
Our numerical results show that the density stratification and magnetic field expansion have 
opposite effects on the oscillating properties of torsional Alfven waves. 
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1 Introduction 



Hannes Alfvcn (1942) predicted the existence of Alfven wave which is one of the magnctohy- 
drodynamic (MHD) waves propagating in magnetized plasmas such as the solar atmosphere. 
Torsional Alfven waves can be observed as temporal and spatial variations in spectral emis- 
sion along the coronal loops (Zaqarashvili 2003). They are an ideal tool for coronal seismology 
as their phase speed depends on plasma quantities within the loop alone, while wave speeds 
of magnetosonic oscillations are influenced by plasma conditions in the ambient medium (Za- 
qarashvili & Murawski 2007). More recently, torsional Alfven waves in the solar atmosphere 
were discovered by Jess et al. (2009) using the high-resolution Swedish Solar Telescope. 

Torsional Alfven modes can have a significant role in coronal heating and solar wind ac- 
celeration, based upon the ability of torsional waves to penetrate easily into the corona (see 
e.g. Ruderman 1999; Copil, Voitenko & Goossens 2008). Possible Alfven wave dissipation pro- 
cesses that may be responsible for heating the solar plasma have been investigated by various 
authors, e.g., dissipation in resonant layers (Pocdts, Goossens & Kcrner 1989, 1990a, b,c; Erdclyi 
& Goossens 1995; Ruderman et al. 1997a,b; Erdelyi 1998; Ruderman 1999; Narain et al. 2001; 
Safari et al. 2006; Karami, Nasiri &; Amiri 2009; Karami &; Bahaxi 2010) and phase mixing (Hey- 
vaerts & Priest 1983; Ruderman, Nakariakov & Roberts 1998; Smith, Tsiklauri Sz Ruderman 
2007; Karami & Ebrahimi 2009). 

People have also paid special attention to the nonlinear effects of torsional Alfven modes. It 
was demonstrated numerically that the observed spiky intensity profiles due to impulsive energy 
releases could be obtained from nonlinear torsional waves (see e.g. Moriyasu et al. 2004; Antolin 
et al. 2008). Taroyan (2009) showed that small-amplitude Alfven waves can be amplified into 
the nonlinear regime by the presence of siphon flows in coronal loops. 

Zaqarashvili Sz Murawski (2007) investigated the evolution of torsional Alfven waves in lon- 
gitudinally inhomogeneous coronal loops. They concluded that the inhomogeneous mass density 
field leads to the reduction of a wave frequency of torsional oscillations, in comparison to that 
estimated from mass density at the loop apex. Also this frequency reduction results from the 
decrease of an average Alfven speed as far as the inhomogeneous loop is denser at its footpoints. 

Copil, Voitenko & Goossens (2010) studied torsional Alfven waves in twisted small scale 
current threads of the solar corona. They showed that the trapped Alfven eigenmodes do exist 
and are localized in thin current threads where the magnetic field is twisted. They pointed out 
that the wave spectrum is discrete in phase velocity, and the number of modes is finite and 
depends on the amount of the magnetic field twist. Also the phase speeds of the modes are 
between the minimum of the Alfven speed in the interior and the exterior Alfven speed. 

Vasheghani Farahani, Nakariakov &: Van Doorsselaere (2010) investigated torsional axisym- 
metric long wavelength MHD modes of solar coronal plasma structures with the use of the second 
order thin flux tube approximation. They concluded that the phase speed of torsional waves 
depends upon the direction of the wave propagation, and also the waves are compressible. 

Verth, Erdelyi &; Goossens (2010) studied the observable properties of torsional Alfven waves 
in both thin and finite-width stratified and expanding magnetic flux tubes. They demonstrated 
that for thin flux tubes, observation of the eigenmodes of torsional Alfven waves can provide 
temperature diagnostics of both the internal and surrounding plasma. They also showed that 
in the finite-width flux tube regime, these waves are the ideal magneto-seismological tool for 
probing radial plasma inhomogeneity in solar waveguides. 

All mentioned in above motivate us to have further investigates on torsional Alfven waves 
by considering the effects of both density stratification and magnetic field expansion on the 
frequencies and eigenfunctions of torsional Alfven modes in the magnetic flux tubes. This paper 
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is organized as follows. In Section 2 we introduce the model and derive the equations of motion. 
In Section 3 we give numerical results. Section 4 is devoted to conclusions. 



2 Model and equations of motion 

We consider an expanding magnetic flux tube of length 2L with longitudinal plasma density as 
typical coronal loop. The tube is assumed to be thin, Va/L <^ 1, where ra is the tube radius 
at the apex. Following Ruderman, Verth & Erdelyi (2008) and Verth & Erdelyi (2008), the 
background magnetic field is assumed to has both radial and axial components with r- and 
z-dependence, i.e. Br = Br{r,z) and B^ = Bz{r,z). The coronal plasma is nearly zero-/^ and 
this yields the magnetic field to be force free. For the selected magnetic field, the electrical 
current is in the (/^-direction. Hence the force free condition, i.e. J x B = 0, is satisfied when 
the electrical current J = V x B = 0. The background magnetic field can be related to a vector 
potential field A as 

B = V X A, (1) 

where 

A=^e,. (2) 

Therefore the radial and axial components of the magnetic field can be expressed in terms of 
the scalar potential V' as 

— Id-tp ldil) 

Br = — Bz = --^-5 (3) 

r oz r or 

where the magnetic field here is perpendicular to 'Vip(r, z), i.e. the magnetic field lines lie in the 
surface 'ip{r, z) = constant. Hence, the equation of the tube boundary is given by ip{r, z) = ipo, 
where Vo is a constant. Here, our aim is to study the torsional Alfvcn waves in which the 
surface of the flux tube tp = ipo has an oscillating motion in the azimuthal direction. If we apply 
the force free condition V x B = to the equilibrium magnetic field (1), we obtain a partial 
differential equation for ijj as 

Verth h Erdelyi (2008) solved the above equation and found the z-component of the magnetic 
field as 



, , cosh(f ) - cosh(l) 



(i-r 

and the radius of the magnetic flux tube boundary as 



cosh(4) — cosh(l 



-1/2 



where B^ ^ = Bz{^L) and = r(ibL) are z-component of the magnetic field and the radius 
of fiux tube at the loop footpoints, respectively. Also F = ^ = f^^jj is the tube expansion 
factor which is defined as ratio of the tube radius at the apex {z = 0) to the tube radius at the 
footpoints {z = ±L). For a tube with constant cross section, the expansion factor is unity but 
for an expanding flux tube we have F > 1. As Ruderman, Verth &; Erdelyi (2008) emphasized 
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the important property of this particular model is that it can describe only magnetic tubes with 
relatively small expansion factors, definitely smaller than 1.87. We also take into account the 
effect of density stratification and assume that the density varies exponentially with the height, 
h, in the atmosphere as p = pfe~H . Here p/ is the density at the footpoints and H is the density 
scale height. Following Ruderman, Verth & Erdelyi (2008) for a half-circle loop with length 2L 
the density can be written as 



p{z) = pf exp 



— 2/x cos 



TTZ 

2L 



L 
'ttH' 



(7) 



where /x is defined as stratification parameter. 

Note that the Alfven velocity in our selected model varies only along the background magnetic 
field. Since we used the thin tube approximation, hence the variation across the magnetic field 
lines is negligible. Ruderman et al. (1997a,b) to study the resonant absorption of torsional 
Alfven waves considered the variation of Alfven velocity only across the magnetic field. 

The linearized MHD equations for a zero-/? plasma are 



Airp 



(V X SB) X B, 



(8) 



(5B = V X X B), 



(9) 



where ^ = (0,0,^,^^) is the Lagrangian displacement of the plasma and (5B = {0,0, SB^) is 
the Eulerian perturbation in the magnetic field. Note that in Eq. (8) due to the force free 
background magnetic field, the term ^^(V x B) x 5B is absent. 
We rewrite Eqs. (8) and (9) in components as 



d%_ 1 
dt^ 47rp 



Br d{r5B^) , ^ d6B^ 

h Bz 



dr 



dz 



(10) 



6B, 



d{BrC<p) , d{B,C^) 



dr 



+ 



dz 



(11) 



Now like Ruderman, Verth & Erdelyi (2008) and Verth, Erdelyi & Goossens (2010) we use a 
non-orthogonal flux coordinate system in which ip becomes an independent variable instead of r, 
i.e. r = r(V', z). In this coordinate system, an arbitrary function /(r, z) will be transformed to 
another function F{ip,z) as f{r,z) = F{'tp{r, z), z). Therefore, using Eq. (3) the r and z partial 
derivatives of / transform to 



(^) =rB,— (^) =^-rBr — 
\dr ' z ^ dip ' ^dzJ'ip dz ^ dip 



(12) 



Differentiating the identities ip = ip{r{ip, z),z^ and r = r{ip{r, z), z^ with respect to z and using 
Eq. (3) one can get 



dr Br 



dr 



dz Bz^ dip rBz 
Using Eqs. (12) and (13) one can rewrite the field Eqs. (10) and (11) as 

d% _ Bz djrSB^) 
dt'^ Airrp dz ' 



(13) 



(14) 
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(15) 



Substituting Eq. (15) into (14) and considering the time-dependence as e the result yields 



where the variables Bz{z), r{z) and p{z) are given by Eqs. (5), (6) and (7), respectively. 
Equation (16) is the same as Eq. (24) in Verth, Erdelyi & Goossens (2010) and also the same as 
the equation for continuum Alfvcn waves on individual magnetic surfaces for a purely poloidal 
magnetic field for any azimuthal wave number m (see Eq. 71 in Poedts, Hermans &; Goossens 



Following Verth, Erdelyi & Goossens (2010) in the thin flux tube limit, all magnetic surfaces 
oscillate with the same frequency, i.e. the frequency is independent of magnetic surface of 
constant as long as ra/L <C 1. Therefore, the torsional Alfven waves described by Eq. (16) 
are global in which the whole loop oscillates with one frequency. For the finite width flux tube, 
Verth, Erdelyi & Goossens (2010) found the local torsional Alfven wave frequencies that vary 
across the magnetic surfaces. 

One notes that the MHD waves have mixed properties. For instance, torsional Alfven waves 
can be coupled to kink waves (see e.g. De Groof Sz Goossens 2000, 2002; Goossens Sz De 
Groof 2001; De Groof, Paes & Goossens 2002; Goossens, Dc Groof & Andries 2002; Goossens, 
Andrics &: Arregui 2006; Goossens 2008). Indeed, the pure torsional Alfven waves require an 
axisymmetric background with a purely poloidal magnetic field (no azimuthal component) and 
axisymmetric motions. 

In the next section we solve the above differential equation using the suitable boundary 
conditions to obtain the eigenvalues uj and the eigenfunctions ^(p{z). 

3 Numerical results 

Here, we solve Eq. (16) using the shooting method to obtain both the eigenfrequencies and 
eigenfunctions of torsional Alfven waves in stratified and expanding magnetic fiux tube. Wc use 
the rigid boundary conditions and assume that ^^{—L) = ^(f){L) = 0. As typical parameters 
for a coronal loop, we assume 2L = 10^ km, pf = 2 x 10~^^ gr cm~^ and B^j = 100 G. 
For such a loop, one finds Alfven speed va = , ^'-^ = 2000 km s"^ and Alfven frequency 

ua:=^ = 0.02 rad s^^ 

The effects of magnetic field expansion and density stratification on the frequencies of the 
fundamental and first-overtone n = 1,2 torsional Alfven modes are displayed in Figs. 1 and 2, 
respectively. Figure 1 shows that for a given stratification parameter /x, the frequencies of the 
fundamental and first-overtone modes decrease with increasing the expansion factor F. This 
result is in agreement with that obtained by Ruderman, Verth & Erdelyi (2008) for the kink 
modes. Figure 2 presents that for a given expansion factor F, the frequencies of the fundamental 
and first-overtone modes increase when the stratification parameter increases. This also is in 
good concord with the result derived by Karami, Nasiri & Amiri (2009) for the kink and fluting 
body waves. One notes that Verth, Erdelyi & Goossens (2010) showed that for a vertical 
stratified and expanding thin magnetic flux tube, for the isothermal case in which the Alfven 
velocity is constant along the magnetic held, the frequency of torsional Alfven waves remains 
unchanged. But for the non-isothermal case in which the Alfven velocity varies along the flux 
tube, the frequency of oscillations for hot and cool tubes decreases and increases, respectively. 




(16) 



1985). 
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The results of Verth, Erdelyi & Goossens (2010) for hot and cool tubes are in good agreement 
with our results presented in Figs. 1 and 2, respectively. In the other words, one can say that 
the cool and hot loops behave like those tubes in which the effect of density stratification and 
magnetic field expansion is dominant, respectively. 

The period ratio P1/P2 of the fundamental and first-overtone n = 1, 2 torsional Alfven modes 
versus the expansion factor and stratification parameter is plotted in Figs. 3 and 4, respectively. 
Note that the period ratio is used as a seismological tool to investigate e.g., longitudinal structure 
(Andries et al. 2005; Andries, Arregui & Goossens 2005) and radial structure (Verth, Erdelyi & 
Goossens 2010) of magnetic loops. Figure 3 reveals that for a given the period ratio P1/P2 
increases with increasing the expansion factor. This is in agreement with the result obtained 
by Verth & Erdelyi (2008) for the kink body modes. Figure 4 illustrates that for a given F, the 
period ratio Pi/ P2 decreases when the stratification parameter increases. This also is in good 
concord with the result reported by Karami, Nasiri & Amiri (2009) for the kink and fluting body 
waves. Note that Verth, Erdelyi & Goossens (2010) pointed out that for a vertical stratified and 
expanding thin magnetic flux tube, the period ratio of torsional Alfven modes for the isothermal 
case is 2. Whereas for the non-isothermal case, the period ratio for both cool and hot tubes 
is bigger than 2. They also considered a semicircular thin magnetic flux tube and showed that 
the period ratio for the isothermal case is 2. But for the non-isothermal case, the period ratio 
for hot and cool tubes becomes bigger and smaller than 2, respectively. The results of Verth, 
Erdelyi &: Goossens (2010) for hot and cool semicircular thin tubes are in good agreement with 
our results illustrated in Figs. 3 and 4, respectively. 

From Figs. 1 to 4 one can conclude that the magnetic field expansion and density stratifica- 
tion have opposite effects on the frequencies and period ratio P1/P2 of torsional Alfven modes. 
This is expected because as a result of magnetic field expansion, Alfven speed becomes smaller 
at the loop apex and the travel time of wave between the two footpoints becomes greater. But 
density stratification causes Alfven speed becomes greater at the loop apex. Hence, the travel 
time of wave between the two footpoints becomes smaller. 

Here, we also investigate the effects of both the magnetic field expansion and density strati- 
fication on the eigenfunctions of torsional Alfven modes. Figures 5 and 6 show the normalized 
anti-node shift Az/L of the first-overtone eigenfunction versus the magnetic expansion factor 
and stratification parameter, respectively. The anti-node shift of the first-overtone mode is one 
of the important quantities in coronal seismology. For instance, it was pointed out by Verth 
(2007) that if the magnetic field becomes weaker towards the loop apex, a useful magneto- 
seismological signature is the anti-node of the 1st harmonic since it shifts towards the loop 
apex. Also Andries, Arregui &: Goossens (2009) used the anti-node shift in the presence of lon- 
gitudinal density stratification as signature of the eigenmode modifications. They showed that 
the anti-node shift of the eigenfunctions of the higher kink-oscillation overtones are even more 
than the first overtone. Figure 5 clears that for a given stratification parameter, the normalized 
anti-node shift increases with increasing the expansion factor. This is in good agreement with 
the result obtained by Verth & Erdelyi (2008) for the kink body modes. Figure 6 presents that 
for a given expansion factor, the normalized anti-node shift decreases when the stratification 
parameter increases. In the other words, the magnetic field expansion causes positive anti-node 
shift and density stratification causes negative anti-node shift. To illustrate this in more detail, 
the eigenfunctions of the first-overtone torsional Alfven modes for the different values of and 
F are displayed in Fig. 7. It shows that as a result of density stratification, the anti-nodes of 
the first-overtone eigenfunction shift away from the loop apex (z = 0). This result was also 
obtained by Andries, Arregui & Goossens (2009) for kink modes. Figure 7 clears also that in 
the presence of magnetic field expansion, the anti-nodes shift towards the loop apex which is in 
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agreement with the result obtained by Verth & Erdelyi (2008) for the kink body modes. It is 
remarkable that Verth, Erdelyi & Goossens (2010) illustrated that for a vertical stratified and 
expanding thin magnetic flux tube, the antinodes of the eigenfunctions for the isothermal case 
remains unshifted. But for the non-isothermal case for a cool tube, the distance between the 
antinodes becomes more spaced out with height and for a hot tube the antinodes become closer 
together with height. The results of Verth, Erdelyi & Goossens (2010) for cool and hot tubes 
are in agreement with our results displayed in the up and down panels of Fig. 7, respectively. 

Note that in some cases in above we compared our results for the frequencies, period ratio 
and eigenfunctions of torsional Alfven modes with those obtained by others for kink modes. 
Although torsional Alfven waves (m = 0) and kink waves (m = 1) have axisymmetric and 
nonaxisymmetric motions, respectively, in the presence of the magnetic field expansion and 
density stratification they show the similar behaviour. This may be caused by this fact that the 
restoring force in both of them is the magnetic tension force. This is why that Goossens et al. 
(2009) used the adjective Alfvenic for kink modes. 

4 Conclusions 

Here, the effects of density stratification and magnetic filed expansion on torsional Alfven waves 
in coronal loops are studied. To do this, a typical coronal loop is considered as an expand- 
ing magnetic flux tube that undergoes a density varying along the tube. The linearized MHD 
equations are reduced to an eigenvalue problem for the azimuthal component of the Lagrangian 
displacement of the plasma. Using the shooting method and under the rigid boundary condi- 
tions for the loop footpoints, both the eigenfrequencies and eigenfunctions of the fundamental 
and first-overtone torsional Alfven modes are obtained. Our numerical results show the following. 

i) For a given density stratification parameter n, the frequencies of the fundamental and 
first-overtone torsional Alfven modes decrease and the period ratio P1/P2 increases when the 
magnetic field expansion factor F increases. 

ii) For a given F, the frequencies of the fundamental and first-overtone torsional Alfven 
modes increase and the period ratio P1/P2 decreases when /x increases. 

iii) Both the density stratification and magnetic field expansion shift the location of the 
anti-nodes of the first-overtone torsional Alfven modes but in the apposite directions. 

All mentioned in above illustrate that the density stratification and magnetic field expansion 
have opposite effects on the oscillating properties of torsional Alfven waves. 
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Torsional Alfven Modes (n=l) 
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Figure 1: Frequencies of the fundamental and first-overtone torsional Alfven modes versus the 
expansion factor, F = for different stratification parameter ^ = (solid line), 2/7r (dashed 

line) and 1 (dash-dotted line). The loop parameters are: 2L = 10^ km, pj = 2 x 10"-^^ gr cm~^, 
Bzj = 100 G. Frequencies are in units of Alfven frequency, ooa = 0.02 rad s~^. 
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Torsional Alfven Modes (n=2) 
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Figure 2: Frequencies of the fundamental and first-overtone torsional Alfven modes versus the 
stratification parameter, ^ = for different expansion factor F = 1 (solid line), 1.35 (dashed 
line) and 1.7 (dash-dotted line). Auxiliary parameters as in Fig. 1. 
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Figure 3: The period ratio P1/P2 of the fundamental and first-overtone torsional Alfven waves 
versus the expansion factor for different stratification parameter fi = (solid line), 2/7r (dashed 
line) and 1 (dash-dotted line). Auxiliary parameters as in Fig. 1. 
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Figure 4: The period ratio P1/P2 of the fundamental and first-overtone torsional Alfven waves 
versus the stratification parameter for different expansion factor T = 1 (solid line), 1.35 (dashed 
line) and 1.7 (dash-dotted line). Auxiliary parameters as in Fig. 1. 
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Figure 5: Normalized shift of the first-overtone anti-node versus the expansion factor for different 
stratification parameter = (solid line) , 2/7r (dashed line) and 1 (dash-dotted line) . Auxiliary- 
parameters as in Fig. 1. 
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Torsional Alfven Modes (n=52) 
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Figure 6: Normalized shift of the first-overtone anti-node versus the stratification parameter 
for different expansion factor T = 1 (soUd hne), 1.35 (dashed hne) and 1.7 (dash-dotted hne). 
Auxihary parameters as in Fig. 1. 
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Figure 7: Eigenfunctions of the first-overtone torsional Alfven modes against fractional length 
C, = z/L for different expansion factor and stratification parameter. In the up panel T = 1.5 
and in the down panel /x = 0.5. Auxiliary parameters as in Fig. 1. 
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